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Low-mass stars, ~ 1 — 2 solar masses, near the Main Sequence are efficient 
at producing 3 He, which they mix into the convective envelope on the giant 
branch and should distribute into the Galaxy by way of envelope loss. This 
process is so efficient that it is difficult to reconcile the low observed cosmic 
abundance of 3 He with the predictions of both stellar and Big Bang nucle- 
osynthesis. In this paper we find, by modeling a red giant with a fully three- 
dimensional hydrodynamic code and a full nucleosynthetic network, that mix- 
ing arises in the supposedly stable and radiative zone between the hydrogen- 
burning shell and the base of the convective envelope. This mixing is due to 
Rayleigh-Taylor instability within a zone just above the hydrogen-burning 
shell, where a nuclear reaction lowers the mean molecular weight slightly. 
Thus we are able to remove the threat that 3 He production in low-mass stars 
poses to the Big Bang nucleosynthesis of 3 He. 
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The standard evolution of a low-mass star (Fig. 1) takes it from a short-lived pre-Main- 
Sequence (PMS) state, in which it contracts and heats up but has not yet become hot enough 
to burn its nuclear fuel, to the long-lived MS state in which slow steady nuclear reactions keep 
the star in thermal equilibrium. After several gigayears (but depending strongly on mass) the 
nuclear fuel is exhausted at and near the center, the star becomes cooler, larger and more lumi- 
nous, and it starts to climb the Red-Giant Branch (RGB). Its outer layers become turbulent and 
convective, and this Surface Convection Zone (SCZ) penetrates deeply into the star; but the SCZ 
is forced to retreat again as the fuel-exhausted core, surrounded by a thin hot nuclear-burning 
shell, advances outwards. During the growing phase the SCZ dredges up and homogenises 
material that, at the earlier MS phase, was processed by nuclear reactions in the interior. 
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Figure 1. Evolution of a low-mass Pop I star in a luminosity-temperature diagram. The model was 
computed in ID, i.e. spherical symmetry was assumed, using the code of (1, 2) with updated equation of 
state, opacity and nuclear reaction rates (3). Surface temperature is in Kelvins, luminosity in Solar units. 

Along the MS stars burn hydrogen in their cores by a combination of the pp chain (in which 
4 protons unite to form a 4 He nucleus) and the CNO tri-cycle (in which the same process is 



2 



catalysed by carbon, nitrogen and oxygen). The former is the more important in low-mass 
stars, 55 1 M , and the latter in more massive stars. However even in the more massive stars 
there is still a shell, somewhat outside the main energy-producing region, where the the pp 
chain partially operates, burning H to 3 He but not beyond. 

Because the pp chain is less sensitive to temperature than the CNO cycle, cores of low- 
mass stars are free of convection; but convective cores develop in higher-mass stars because 
CNO energy production is too temperature sensitive for radiation to stably transport the energy. 
Above ~ 2 M this convective core is large enough that the 3 He produced is convected into the 
center of the star and burnt there. But in stars of lower mass 3 He accumulates (4) in a broad 
zone outside the main energy-producing region (Fig. 2a). 3 He is enriched above its assumed 
initial value (2.10 -4 by mass (5), the same as its surface value in this plot) in a broad peak 
extending over nearly half the mass (as well as about half the radius) of the star. The maximum 
3 He abundance in this peak is a factor of ~ 18 larger than the initial value. 
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Figure 2. (a) Profiles of the abundances of certain nuclei in a star which has evolved to roughly the 
end of the MS (Fig. 1; T~ 5000 K, L ~ 2 L ). X H is orange, 4 He is red, 16 is blue, 12 C is black, 14 N 
is green and 3 He is yellow. 3 He shows a major peak where the abundance reaches ~ 18 times the initial 
(surface) abundance, (b) The same star later, when the SCZ reaches its maximum inward extent (Fig 1). 
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The 3 He peak has been homogenized, to a factor of 8 larger than its initial value. The inert, H-depleted 
core is about 0.27 M & . 

On the lower part of the RGB (Fig. 1), a large SCZ develops, which mixes and homogenises 
the outer ~ 0.7 M (Fig. 2b). The surface 3 He abundance is raised from the initial of 2.10 -4 to 
~ 1.6.10~ 3 , i.e. by a factor of ~ 8. 

As the star climbs the RGB beyond the point (Fig. 1) where the SCZ penetrates most deeply, 
the SCZ is diminished by (a) nuclear burning below its base, in a zone that marches outwards, 
and (b) stellar-wind mass loss from its surface. The evidence for the latter is that the next long- 
lived stage after the RGB is the Horizontal Branch (HB), and HB stars appear to have masses 
that are typically 0.5 — 0.6 M Q , substantially less than the masses of stars capable of evolving 
to the RGB in less than a Hubble time (6, 7). Process (b) leads to enrichment of the interstellar 
medium (ISM) in 3 He (8, 9, 10). 

Yet the ISM's abundance of 3 He, at ~ 5.10~ 5 by mass, is little different from that predicted 
by Big Bang nucleosynthesis. This is a major problem (11, 12): either the Big Bang value is 
too high, or the evolution of low-mass stars is wrong. 

In this paper we identify a mechanism by which low-mass stars destroy (on the RGB) the 
3 He that they produced during their MS evolution. Although we illustrate this with a star like 
the Sun, regarding both mass and initial composition, we emphasise that exactly the same ap- 
plies to low-mass metal-poor stars ('Population II'), which may have been more important than 
metal-rich ('Population F) stars like the Sun throughout the earlier part of Galactic history in 
determining the 3 He abundance of the interstellar medium. The process is largely independent 
of mass provided it is fairly low: 1—2 M Q for Pop I and 0.8 — 1.6 M & for Pop II. 

Once the SCZ has reached its deepest extent, part- way up from the base of the RGB, it re- 
treats, and can be expected to leave behind a region of uniform composition with 3 He enhanced 
(Fig. 2b). This region is stable to convection according to the usual criterion that the tempera- 
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ture gradient should be sub-adiabatic, and is quite extensive in radius although small in mass. 
The H-burning front moves outwards into the stable region, but preceding the H-burning region 
proper is a narrow region, usually thought unimportant, in which the 3 He burns. The reaction 
that mainly consumes it is 3 He ( 3 He, 2p) 4 He , which is an unusual reaction in stellar terms be- 
cause it lowers the mean molecular weight: two nuclei become three nuclei, and the mean mass 
per nucleus decreases from 3 to 2. Since the molecular weight (//) is the mean mass per nucleus, 
but including also the much larger abundances of X H and 4 He that are already there and not tak- 
ing part in this reaction, this leads to a small inversion in the /x-gradient. The inversion is tiny 
(Fig. 3): it is in about the fourth decimal place. But our 3D modeling shows the inversion to be 
hydrodynamically unstable, as we should expect from the classic Rayleigh-Taylor instability. 
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Figure 3. The profile of reciprocal molecular weight (1///) as a function of mass in solar units, at 
three successive times (red, then green 2 Myr later, then blue 2Myr later still). 
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At a stage (Fig. 3) when the SCZ has only just begun to retreat, there is no bump in 1/fj,, but 
just a slight distortion at about 0.248 M . This is because the 3 He consumption is taking place 
in a region where there is still a substantial /i gradient left over from earlier history. But as the 
H-burning shell moves out (in mass), the 3 He-burning shell preceding it moves into a region of 
more uniform 1 H/ 4 He ratio, and so the peak in 1/ ji begins to stand out. By the time the leading 
edge of the shell has moved to 0.25 M there is a clear local maximum in which persists 
indefinitely as the H-burning shell advances and the convective envelope retreats. 

At a point somewhat beyond this in the evolution of our ID star (Fig. 1) we mapped the ID 
model on to a 3D model and used the hydrodynamic code 'Djehuty' developed at the Lawrence 
Livermore National Laboratory (13, 14, 15). The code is described most fully in the third of 
these papers. Although Djehuty is designed to deal with an entire star, from center to pho- 
tosphere, we economised on meshpoints by considering only the region below the SCZ. It is 
important for numerical purposes that the ID and 3D codes use exactly the same approxima- 
tions for physical processes, e.g. equation of state, nuclear reaction rates, opacities. 

The location of the starting model of the 3D calculation is shown on Fig. 1. If we had been 
clear before starting the 3D calculation that the 1/ fi bump was going to cause mixing, we would 
have started further down, at the point where the bump first presents itself, which is just above 
the point labelled 'deepest penetration'. It has become clear that our unexpected mixing will 
begin around here, and in practice we expect (as discussed below) that almost all of the 3 He 
in the SCZ will have been consumed by the time the model reaches the point where our 3D 
calculation started. Since 3D modeling is very expensive of computer time, we have chosen 
not to redo the calculation for an earlier starting-point. Fig. 4 is a cross-section of the starting 
model for the 3D run, and shows the yU-inversion as a ring well outside the burning shell. 
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Peak in 1/M-from 3 He burning 




Figure 4. A color-coded plot of ^ on a cross-section through the initial 3D model. The shell where 
the /Li-inversion occurs is the yellow region sandwiched between a yellow-green and a rather darker green. 
The inversion is at a radius of ~ 5.10 7 m. The base of the SCZ is at ~ 2.10 9 m, well outside the frame, 
and the surface of the star is at ~ 2.10 10 m. . 

Following the early development of the initially-spherical shell on which has a con- 
stant value near its peak (Fig. 5), the surface has begun to dimple after only ~ 800 sees, and 
by 2118 sees the dimpling is very marked and the surface has begun to tear. Some points have 
moved ~2% radially, i.e. ~ 10 6 m, indicating velocities of ~500m/s. The mean velocity de- 
creases slightly in the passage from the second to the fourth panel. Other spherical shells, well 
away from the inversion on either side, show no such dimpling, at least until the influence of 
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the inversion has spread to them. A movie of which Fig. 5 is four frames is given as Movie. SI 
in Supplementary Online Material. 




Figure 5. The development with time of a contour surface of mean molecular weight near the peak in 
the blue curve of Fig. 3. The contour dimples, and begins to break up, on a timescale of only ~ 2000 sec. 

The velocity we see is roughly consistent with the expectation that it should be v 2 ~ glAfj,/ //, 
where g is the local gravity and / is the local pressure scale height. The motion appears turbulent, 
and has the effect of diluting the inverse molecular-weight gradient, but it cannot eliminate it. As 
the turbulent region entrains more of the normally stable region outside it yet below the normal 
convective envelope, it brings in fresh 3 He, which burns at the base of this mixing region, thus 
sustaining the inverse molecular- weight gradient. Ultimately this turbulent region will extend 
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to unite with the normally-convective envelope, so that the considerable reservoir of 3 He there 
will also be depleted. If its speed of ~ 500 m/s is maintained the time for processed material to 
reach the classically unstable SCZ is only about 1 month, while the time for the H-burning shell 
to burn through the ~ 0.02 M Q layer is over 10 6 yrs. 

The above argument establishes that the mixing in the SCZ is extended below the classi- 
cal convective limit, and that it is very fast compared to the nuclear timescales of either the 
hydrogen-burning shell or the 3 He-burning reaction. We estimate from the nuclear-burning 
rates we find that as the hydrogen shell burns outwards the 3 He will be destroyed in 16 times as 
much mass as the hydrogen shell burns through. 

We believe that the extra mixing that we have seen gives a satisfactory answer to the prob- 
lem of matching the 3 He abundance of Big Bang nucleosynthesis. Although low-mass stars 
do indeed produce considerable amounts of 3 He on the MS, this will all be destroyed by the 
substantially deeper mixing that we now expect on the RGB. 

Our deeper mixing can also be relevant to further problems. According to the classical 
models of RGB stars, there is no further modification to the composition in an RGB convective 
envelope after it has reached its maximum extent early on the RGB. Yet observations persis- 
tently suggest that the ratios 13 C/ 12 C and 14 N/ 12 C both increase appreciably as one goes up the 
RGB (16, 17). Both these ratios can be expected to increase only if the material in the envelope 
is somehow being processed near the H-burning shell. Our model makes this very likely. Al- 
though the yu-inversion that we find is somewhat above the main part of the H-burning shell, it 
is not far above and we can expect some modest processing of 12 C to 13 C and 14 N. According 
to (18), it appears to be necessary for some extra mixing to take place beyond the point on the 
RGB where the SCZ has penetrated most deeply; that is exactly the point where our mechanism 
should start to operate. In (18, 19, 20, 21) it was suggested that rotation in the region between 
the SCZ and the hydrogen-burning shell might be responsible for the required mixing. We do 
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not dispute the possible importance of rotation; however we emphasise that the mechanism we 
have discovered is not ad hoc, but simply arises naturally when the modeling is done in 3D. 
This mixing occurs regardless of possible variables like rotation and magnetic fields. It seems 
possible to us that different rates of rotation might vary the efficiency of our process, and we 
intend to investigate models with rotation in the future. 

Correlations between abundance excesses and deficits of various elements and isotopes in 
the low-mass evolved stars of globular clusters have been reported in (16). Although it is hard 
to distinguish star-to-star variations that may be due to evolution from those that may be due to 
primordial variation, we expect our mechanism to lead to substantial evolutionary variations. 

We feel that our investigation demonstrates particularly clearly the virtue of attempting to 
model in 3D, where the motion evolved naturally, and to a magnitude that initially surprised us. 
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Figure 1. Evolution of a low-mass Pop I star in a luminosity-temperature diagram. The model was 
computed in ID, i.e. spherical symmetry was assumed, using the code of (1, 2) with updated equation of 
state, opacity and nuclear reaction rates (3). Surface temperature is in Kelvins, luminosity in Solar units. 

Figure 2. (a) Profiles of the abundances of certain nuclei in a star which has evolved to roughly the 
end of the MS (Fig. 1; T ~ 5000 K, L ~ 2 L Q ). X H is orange, 4 He is red, 16 is blue, 12 C is black, 14 N 
is green and 3 He is yellow. 3 He shows a major peak where the abundance reaches ~ 18 times the initial 
(surface) abundance, (b) The same star later, when the SCZ reaches its maximum inward extent (Fig 1). 
The 3 He peak has been homogenized, to a factor of 8 larger than its initial value. The inert, H-depleted 
core is about 0.27 M Q . 

Figure 3. The profile of reciprocal molecular weight as a function of mass in solar units, at 

three successive times (red, then green 2 Myr later, then blue 2Myr later still). 

Figure 4. A color-coded plot of fi on a cross-section through the initial 3D model. The shell where 
the /n-inversion occurs is the yellow region sandwiched between a yellow-green and a rather darker green. 
The inversion is at a radius of ~ 5.10 7 m. The base of the SCZ is at ~ 2.10 9 m, well outside the frame, 
and the surface of the star is at ~ 2.10 10 m. . 

Figure 5. The development with time of a contour surface of mean molecular weight near the peak in 
the blue curve of Fig. 3. The contour dimples, and begins to break up, on a timescale of only ~ 2000 sec. 
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